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Even for the observed luminosity distance Dl{z), which suggests the existence of dark 
energy, we show that an inhomogeneous dust universe solution without dark energy is possible 
in general. Future observation of Dl{z) for 1 < 2 < 1.7 may conhrm or refute this possibility. 


§1. Introduction 

Recent measurements of the luminosity distance Dl(z) using Type la supernovae 
suggest that an accurate value of Dl{z) may be obtained in the near future. 
In particular, SNAP^^ should provide us the luminosity distances of ~2000 Type la 
super novae with an accuracy of a few percent up to z ~ 1.7 every year. Also, from 
observation of the first Doppler peak of the anisotropy of the CMB, it is suggested 
that the universe is flat,®^’®^ and this may be proved in the future from obervations 
by MAP and Planck. Under the assumption of the homogeneity and isotropy of 
our universe, these observations suggest that dark energy is dominant at present. 
In an attempt to determine the nature of dark energy, many arguments have been 
given. Recently, some mechanisms to account for the observed tiny but finite dark 
energy are proposed. However, at present we do not have a firm and reliable 
theoretical basis to investigate such a small energy scale compared with the Planck 
scale. In short, the nature of dark energy under the assumption of the homogeneity 
and isotropy of our universe is still a great mystery. 

From the observed isotropy of the CMB, assuming that we are not in a special 
part of our universe, the universe should be homogeneous. However, if our posi¬ 
tion in the universe is special, the universe might be inhomogeneous, although the 
CMB is isotropic. Such cosmological models have been constructed using spherically 
symmetric models in which we are near the symmetric center. Some authors have 
considered such models to interpret the SNIa data for small z, as well as large 
z assuming a void structure to avoid dark energy. Such possibilities may be 

regarded as absurd. However, our point of view in this paper is to construct a possi¬ 
ble inhomogeneous dust universe derived from the observed Dl(z). If such a model 
is consistent with present observational results, the inhomogeneous universe should 
be examined more seriously, because the dark energy solution is also absurd in the 
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sense that it is ~ 120 orders of magnitude smaller than the Planck scale. In short, 
we suppose that the question we need to answer to roughly reduces to the following: 
Which is more absurd, dark energy or an inhomogeneous universe? In the former 
case, there is no reliable theory to examine the problem at present, while the latter 
case can be studied in the frame-work of known theories. We would like to point 
out that it is not taste but, rather, future observations that will confirm either dark 
energy or an inhomogeneous universe. 

The analysis of high redshift supernovae gives us the luminosity distance-redshift 
relation Dl{z) along the observational past null cone up to z ~ 1.1) 3) data 
fit well with Dl{z) in the homogeneous and isotropic universe with 17^ = 0.3 and 
Qa = 0.7 given by 


Dl{z) 


J_. . N __ 

Hq Jo ^ j7m(l + + [2a 


(M) 


In this paper we assume that Dl{z) is given by Eq. (IT) with 17^ = 0.3 and 
17/1 = 0.7 for z ^ 1. This is done for the sake of simplicity to make the arguments 
clearer. In particular, we do not wish to claim that Dl{z) with f2m = 0.3 and 
17/1 = 0.7 has been confirmed. While Dl{z) for 1 < z < 1.7 is not certain even at 
present and will be obtained in the future, for example, by SNAP. Since the scale 
factor a obeys 

d dvr, , , , 

- = -^(p + 3p, 1-2 

a 6 

Di{z) with Qrn = 0.3 and l7/i = 0.7 implies that the present universe is accelerating, 
while for the dust universe {p = 0), a should be decelerating. Therefore it may 
be concluded that observations is inconsistent with the inhomogeneous dust model. 
However, the point is that to determine Dl{z), we are observing Type la supernova 
events that occurred at past times in spatial positions separated from us. In the 
inhomogeneous universe model, the time dependence of a at a point separated from 
us differs from that of our position, so that we may obtain an apparent accelerating 
universe even though the dust universe is decelerating locally. 

Before ending this introduction, we comment on some other works relevant to 
this paper. The inhomogeneous scenario is not the only alternative to dark energy. 
Giving up the assumption that the cosmic substratum is composed of perfect fluids, 
bulk pressures that differ from the kinetic pressure can be allowed. The assumption 
of an effective anti-friction force leads to a model that has only one dark component 
(CDM) and is consistent with the CMB and SNIa data. Also, there is an approach 
somewhat related to that presented in this paper (though with different motivation) 
that has been used. 


§2. Formulation 


The line element of a spherically symmetric dust universe is given by 

{R'{t,r)r 


ds = —dt + 


-dr^ + R^{t, r)dQ^, 


1 -|- 2E{r)r‘^ 


(2T) 
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where the prime indicates differentiation with respect to r. The solution to the 
Einstein equations is known as the Lemaitre-Tolman-Bondi (LTB) spacetime, 
given by 


R = 


4:7rp{t, r) = 


2M(r) 

M' 

Wr'' 


+ 2E{r)r‘^, 


( 2 - 2 ) 

(2-3) 


where the dot indicates differentiation with respect to t. The solution of Eq. (2-2) is 
given by 

M M 

R{t,r)= (f{t,r), t-tB{r) =^{t,r)- 


where 


and 


<P = 


e{r)r‘^ 


err = 


cosh ?7 — 1, 
2 ’ 


(e(r)r^)^/^ 


2E(r)r^, (-E(r) > 0) 

1, (E(r) = 0) 

—2E{r)r‘^, {E{r) < 0) 


sinh?] — 7], {E{r) > 0) 

« = i f, (E(r) = 0) 

rj — sinr/. (E'(r) < 0) 


(2-4) 


(2-5) 


( 2 - 6 ) 


1 — cos 7], 

In the general solutions of the LTB models, there are three arbitrary functions 
M{r),E{r) and tsir). M{r) is regarded as the gravitational mass function, and 
we can set M{r) = Mgr^, redefining r. tsir) corresponds to the local BigBang time. 
E(r) determines the local curvature radius or the local specific energy. The functions 
tsir) and E{r) should be chosen to reproduce the observed Dl{z). This means that 
we have only one constraint for two arbitrary functions. 

The observational past null cone is specified in the form, t = t{r). We denote 
the areal radius R out = t{r) by TZ. Then, by Eq. (2-2), we can regard R out = t{r) 
as a function oilZ, E and r: 


i?(t»,r) =7^l(7^,E,r) = J^^^ + 2Er^. 


(2-7) 


By differentiating Eqs. (2-2) and (2-4), R' and R' on t = t{r) can be expressed as 
functions of TZ, i, E, E', ts, t'^ and r:*) 


R'{i{r),r) = 7^2 (tZ, i, E, E',tB,tB,r"j 


= -{TZ--, 
V 2 L 

TZ 

4—) 

r 


t — tB 


^ “ RR'b 


If E{r) = 0, we should omit the terms proportional to E' in Eqs. (2-8) and (2-9). 


( 2 - 8 ) 
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and 


R'{i{r),r) = 7^3 (jZ, i, E, E'NbNbn) 


1 / Mor^ , 

Mor^ 7^l 

+ + —• 

The observational past null cone t = t{r) satisfies 

di _ Tl2{E,i,E,E'N bN'bN) 

dr \/l + 2Er‘^ 

The redshift z{r) along the past null cone is given by 
dz 1 + z 
dr Vl + 2£'r2 ' 

The total derivative of TZ on the past null cone is written 


E' 

'e 


=7^3(7^, i, E, E', tB,t'B, r) 


dTZ 7Zi(7Z, E,r) 

dr 


Vl + 2Er^ 


TZ2{'JZ,i,E,E',tB,t'BN)- 


(2-9) 


( 2 - 10 ) 


( 2 - 11 ) 


( 2 - 12 ) 


Our basic equations are Eqs. (2T0) —(2T2). These three equations can be regarded 
as a system of first-order ordinary differential equations for three of the five functions 
7Z{r), t{r), E{r), tsir) and z{r). In order to integrate these equations, we must 
specify two conditions on these five functions. The luminosity distance Dl{z) is 
related to 1Z‘^^'> as 

7 ^ _ Rl{z) (2-13) 

(1+Z)2- ^ 

As mentioned above, we assume that Dl{z) is given by Eq. (IT). We will specify 
one further condition on E, or a combination of them. 


§3. Results 


3.1. Results of the BigBang time inhomogeneity 

We first consider a pure BigBang time inhomogeneity. In this case, the curvature 
function E{r) is set to a constant value. Erom Eqs. (2T1) and (2T2), we have 
equations for z{r) and the Big-Bang time function tsir). The model is specified by 
l7o = ‘^Mq/Hq, which is the present central density, SMq/Att, divided by the present 
central critical density, pcrit = 3Rg/87r, where Hq is the present central Hubble 
parameter, and we set it to unity. We numerically integrated these two differential 
equations from r = 0 for ten l7o from 0.1 to 1.0. The initial conditions are given by 
z = 0 and ts = 0- 

From Eq. (2-3), R' > 0 for positive density, while from Eq. (2T1), R' > 0 for 
monotonically increasing z{r), so that the integration is terminated when either of 
the inequalities 


R' > 0 or R' > 0 


(3T) 



Is dark energy the only solution to the apparent aceeleration of the present universe^ 

is violated. In Fig. 1, we display the relation between the parameter l7o and the 
redshift at the time that the integration is terminated. For small value of Hq, i7o = 
0.1 - 0.4 (open triangles), shell-crossing singularities appear when dlZ/dz = 0. For 
large value of l7o, l?o = 0.5 - 1.0 (open square), the second condition in Eq. (3-1) is 
violated first. This occurs when TZ = 2M. 
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<C 
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Fig. 1. Plots of the maximum redshift of the time that either of the inequalities in Eq. (3-1) 
is violated as a function of the present density parameter. The open triangles and the open 
squares correspond to the BigBang time inhomogeneity. The crosses correspond to the curvature 
inhomogeneity. 
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Figure 2 plots the BigBang time functions ts for each Qq. For all Hq, the 
BigBang time functions ts decrease as z increases up to z ~ 0.5. This result is 
related to the fact that the expansion of our universe appears to be accelerating up 
to z ~ 0.5. In inhomogeneous models, an apparent acceleration is realized if the 
recession velocity of mass shells does not increase rapidly along the observational 
past null cone as in the case of a homogeneous and isotropic universe filled with 
dust. To construct such a situation in our model, we need to prepare an older shell, 
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i.e., one that is more decelerated by gravity, for more distant shells on the past null 
cone. This is the reason that the function ts decreases. 

In Fig. 3 we plot the redshift space density, 


AnR^R'dr 


7-2 dj. 

^^0 ~~n ~i~ PcYlt 

dz 


(3-2) 


along the past null cone. Observations of the mass distribution along the past null 
cone would give us this density profile. 



Fig. 3. Plots of the redshift space density p divided by the central critical density. The dotted 
curve represents the flm = 0.3, Qa = 0.7 homogeneous model. 

3.2. Results of curvature inhomogeneity 

Next, we consider the pure curvature inhomogeneity. In this case the BigBang 
time function t_B(r) is set to zero. From Eqs. (2T0) ~ (2T2) we obtain three differen¬ 
tial equations for the three variables z{r), E{r) and t{r). We numerically integrated 
these three differential equations from r = 0. The initial conditions are given by 
z = t), E = {1 — f7o)/2 and 


where 


t(0) 


f7o (sinhT/o - r/p) 

2 (l-f2o)^ 



(3-3) 


(3-4) 


The present central cosmological time t(0) and ryo are obtained from Eq. (2-4). 

In Eig. 1, we show the relation between the parameter Qq and the redshift at 
the time that the integration is terminated (cross marks). Eor the case of a curvature 
inhomogeneity, it is found that the second condition in Eq. (3T) is violated first. 

Eigure 4 displays the curvature functions E for variations values of Qq. We can 
see E decreases as z increases, except in the i?o = 1-0 case. 
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The result that E decreases is consistent with the apparent acceleration. The 
value of E determines the specific energy of the dust elements, so that the “initial” 
velocity is smaller for more distant shells. This causes apparent acceleration, because 
the velocity at r = 0 can be largest. 


0.4 


0.2 
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Fig. 4. Plots of the curvature functions. 

Figure 5 displays the redshift space density p along the past null cone as a 
function of z. 
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I 
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Fig. 5. Plots of the redshift space density divided by the present central critical density. The dotted 
curve represents the flm = 0.3, Qa = 0.7 homogeneous model. 

§4. Summary and discussion 

In this paper we have constructed inhomogeneous dust models without dark 
energy. We find that these models are consistent with the observed Dl{z) up to z = 
1, as from Fig. 1 no difficulties are encountered up to z ~ 1 for any set of parameter 
values in both the BigBang time inhomogeneity and curvature inhomogeneity cases. 
For z > 1, we have difficulties in our inhomogeneous dust models. Recently, the SNIa 
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at a redshift of ~ 1.7 was found with rather large uncertainties. However, 

only a single SNIa at a redshift of ~ 1.7 is not enough to construct an accurate 
Dl{z), although that result seems to rule out the ‘grey-dust’ hypothesis. In addition, 
the results of a recent investigation of the effect of gravitational lensing on this 
SNIa suggests that the grey-dust model may be consistent with the observational 
data. If future observations confirm Dl{z) up to z ~ 2 with 17^ ~ 0.3 and 

~ 0.7 , it can be concluded that our inhomogeneous dust models are incompatible 
with the observations and that some form of dark energy is likely to exist. However, if 
future observations confirm that Dl{z) for z > 1 is not consistent with Eq. (1-1), the 
plausibility of our inhomogeneous dust models should be studied more extensively. 
In such a case, the first Doppler peak as well as the higher ones will give us another 
constraints on the inhomogeneous universe models. 

It may believed that the existing observations for 0 < z < 1, such as (i) evolution 
of cluster abundance, (ii) lensing rate, and (iii) ages of stellar populations, already 
rule out the inhomogeneous models. 

Using the cluster temperature evolution data for 0.3 < z < 0.8, it was found 
that the best-fit value for is Om = 0.45 ±0.1 for open universe and Qm = 
0.3±0.1 for flat universe. However, recent analysis shows that the systematic error 
is comparable to the statistical error. Therefore, we may say that 0.1 < Qm < 0.5 
for 0.3 < z < 0.8 data. It is not clear whether the Press-Schechter formalism can be 
applied to our inhomogeneous models. One possible estimate of cluster abundance 
could be based on the locally homogeneous approximation. As we know, massive 
cluster evolution is very sensitive to matter density. It seems that a model whose 
local density parameter differs greatly from the best-fit value would not be able to 
explain the observed cluster evolution. The pure curvature inhomogeneity case with 
f7o ~ 0.2 may not survive, because it is approximated by a flat universe at large z. 
Also the BigBang time inhomogeneity case with f7o ~ 1-0 cannot survive. However, 
it can be expected that the pure BigBang time inhomogeneity with Qq ~ 0.5 and 
the pure curvature inhomogeneity with Qq ~ 0.1 will predict the observed cluster 
abundances. 

The estimate of the lensing rate and the distribution of the separation of the 
images depend on the model used for the mass distribution of the lensing object and 
the luminosity function of the source objects as well as the cosmological parameters. 
However, it has been shown that the dependence on the lens model and parameters 
is much stronger than that on the cosmological parameters. In addition, the 

mass distribution of the lensing objects should depend strongly on the baryon density 
Therefore, we conclude that the estimate of the cosmological parameters from 
the lensing rate and the distribution of the separation of the images is difficult at 
present, and for this reason, we cannot rule out the inhomogeneous model. 

As shown in Fig. 6, the look back times along the past null cone differ little 
between the inhomogeneous model and the corresponding homogeneous model with 
cosmological constant for z < 0.5. For z ~ 1, a difference appears, but some of 
the inhomogeneous models do not differ greatly from the homogeneous model even 
in that case. The ages of a stellar population could not be used to distinguish the 
inhomogeneous model from the homogeneous one. 
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Fig. 6. Plots of the look back time along the past null cone. The solid curve represents the ho¬ 
mogeneous fim = 0.3, I2a = 0.7 case. The broken and dotted lines denote the pure BigBang 
time and the pure curvature inhomogeneity cases with 17o = 0.1, 0.3, 0.7 in descending order, 
respectively. 


We have found that the model dependence, including various undetermined pa¬ 
rameters, and the observational uncertainty are much larger than the dependence on 
the cosmological parameters. Therefore we believe that these observations cannot 
easily rule out the inhomogeneous model. 

Before finishing, we give a brief comment on how we can be positioned away 
from the center of the symmetry. A displacement from the center would correspond 
to a dipole mode of CMB. Therefore we can be positioned ~ 50 Mpc away from the 
center. 

In conclusion, dark energy is not the only possible solution of the apparent 
acceleration of the present universe, as inhomogeneous dust models can also account 
for current observations. 
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